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Abstract. The paper deals with the existence and almost periodic homog- 
enization of some model of generalized Navier-Stokes equations. We first es- 
tablish an existence result for non-stationary Ladyzhcnskaya equations with 
a given non constant density and an external force depending nonlinearly on 
the velocity. Next, the density of the fluid being non constant, we combine 
some compactness arguments with the sigma-convergence method to study the 
asymptotic behavior of the velocity field. 



1. Introduction 

The Navier-Stokes equations model the motion of Newtonian fluids. In order to 
understand the phenomenon of turbulence related to the motion of a fluid, several 
mathematical models have been developed and studied over the years. We refer to, 
e.g. [1 EH H ES3 HU, just to cite a few. 

To investigate the turbulence in non-Newtonian fluid, we consider a model close 
to the Ladyzhenskaya one [20] . In this model, the viscosity and the density are non 
constant. To be more precise, let e > be a small parameter representing the scale 
of the inhomogeneities. The equation of the motion reads 
(f.I) 

P e -^T - div (a E Vu e + b £ |Vm £ | p ~ 2 V« e ) + (u e • V)m £ + Vg £ = p £ f e {-,u e ) in Q T 

where p is the non constant positive known density which is bounded from above 
and from below away from zero, u £ and q e are the unknown velocity and pressure, 
respectively. The viscosities a = {aij)\<i.j<N arid b both depend on the spatial and 
time variables. The matrix a is coercive whereas the function b is positive, bounded 
from above and from below away from zero. The density and the viscosities are 
scaled as follows: p s {x) = p (| ) , a £ Jx,t) = (§> p) and6 e (x,i) = &(f,p) 
whereas the external force f £ (-,u £ ) which is a Lipschitz function of the velocity u £ 
is scaled as follows r)(t) = f (4r>J") for (x, t) £ Qt and r £ R N . The problem 
is stated in details in the beginning of section 2. 

The equation (jl.ll) models various types of motion of non-Newtonian fluids. We 
cite a few examples. If a = (2A)#y)i<i,j<jv (<% the Kronecker delta) and b = V\ 
where vq and v\ are positive constants, then we get the Ladyzhenskaya equations. 
In that case, the analysis conducted in [IS] reveals that one may either let v\ — >• 
(and get the usual Navier-Stokes equations) or let v — > 0, and hence we are led to 
the power-law fluids equations. In particular when p = 3, we get the Smagorinski's 
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model of turbulence |36] with vq being the molecular viscosity and v\ the turbulent 
viscosity. Another model included in (jl.ljl is the equation of incompressible bipolar 
fluids 0. 

There are several works dealing with the Navier- Stokes equations with external 
force depending on the velocity. These types of equations are commonly known as 
the generalized Navier-Stokes equations. This work is different from the previous 
ones in the sense that it combines both a non-constant density, the non-Newtonian 
fluid effect, and the external force depending on the velocity. It should be noted that 
if in equation (11.11) we replace the gradient Vit by its symmetric part |(Vm + V t u), 
then thanks to the Korn's inequality, the mathematical analysis does not change, 
although the model becomes in this case, physical. 

Started in the 70 's the mathematical theory of Homogenization is nowadays 
divided into two major components: the individual homogenization theory (also 
known as the deterministic homogenization theory) and the random homogeniza- 
tion theory (also known as the stochastic homogenization theory). In this paper we 
are concerned with a special case of individual homogenization theory, namely the 
almost periodic one. More precisely, we assume throughout the paper that the den- 
sity p, the viscosities a and b and the source term / are almost periodic functions as 
specified in the beginning of Section 4. It should be noted that we could consider a 
more general deterministic setting with coefficients satisfying general deterministic 
assumptions covering a large set of concrete behaviors (see e.g., [131 H31 [53]) such 
as the periodic one, the almost periodic one, the convergence at infinity and many 
more besides. We only deal with almost periodic homogenization for the sake of 
simplicity. Our results carry over mutatis mutandis to the general deterministic 
setting. Concerning the homogenization of (II. ip . it should also be stressed that 
this is the first time that such an analysis is conducted beyond the periodic setting. 

The goal of this paper is twofold. We first establish an existence result for 
(|l.ip and then perform the homogenization process for (|l.ip in the almost periodic 
setting. The space dimension is either 2 or 3. The only place where the analysis 
changes depending on the dimension is the proof of the inequality (|2.8[) which is very 
useful in the proofs of the existence result, theorem [T] and the compactness result, 
Proposition [TJ We stress that the existence result Theorem [T] is not solving the 
millennium problem of the existence and smoothness of the Navier-Stokes equations 
as stated by Fefferman |17j since we only consider weak solutions here with no 
regularity result. 

The paper is organized as follows. Section 2 deals with the complete statement of 
the problem, the proof of the existence result and some a priori estimates. In Section 
3 we gather some necessary tools about the ^-convergence method (which is just 
the appropriate generalization of the well-known two-scale convergence method) in 
the algebra of continuous almost periodic functions. Finally in Section 4, we state 
and prove the homogenization result. 

Throughout Section 3, vector spaces are assumed to be complex vector spaces, 
and scalar functions are assumed to take complex values. We shall always assume 
that the numerical space M. m (integer m > 1) and its open sets are each equipped 
with the Lebesgue measure dx = dx\...dx m - 



HOMOGENIZATION OF LADYZHENSKAYA EQUATIONS 



3 



2. Statement of the problem: existence result and a priori estimates 

2.1. Problem setting: existence result. We consider A-dimensional problem, 
N = 2,3. In what follows, all the function spaces are real- valued spaces and scalar 
functions assume real values. 

Let 1 + < p < oo, and T > a real number. Let Qt = Q x (0,T) where 
Q is a bounded smooth domain in R . We consider the following well-known 
spaces [221 EH]: V = {(p £ C$°(Q) N : div ip = 0}; V = closure of V in VF 1 ^(Q) Ar ; 
H = closure of V in L 2 (Q) N . In view of the smoothness of Q, it is known that 

V = {u e Wq'^Q)^ : divw = 0} and H = {it e L 2 (Q) N : divw = and 
u| a g -ti = 0}, where denotes the trace of it on and n is the outward unit 
vector normal to dQ. The space V is hence endowed with the gradient norm and 
H with the L 2 -norm. For the sake of simplicity, we denote the respective norms of 

V and H by ||-|| and |-|. We denote by (, ) the inner product in H, as well as the 
scalar product in Mr. The associated Euclidean norm in M. N is also denoted by |-|. 
All duality pairing will be denoted by (, ) without any specification of the spaces 
involved. 

With all this in mind, we shall consider the functions a, 6, p and / constrained 
as follows. 

(Al) The matrix a = (fflij)i<i,j<JV G L oc (R N+1 ) NxN satisfies = a 3i and 

(a(y, t)A) • A > v |A| 2 for all Ael" and a.e. (y, r) e R N+1 

for some positive vq. 
(A2) The function b £ L 00 (IR Ar+1 ) is such that v\ < b < v 2 for some positive v\ 
and V2- 

(A3) The density function p belongs to e L°°(R N ) and satisfies A" 1 < p(y) < A 

a.e. y £ R N , for some positive A. 
(A4) The mapping / : 1 x -> R N , (t,r) h4 f(t,r) is continuous with the 

following properties: 

(i) / maps continuously Kxl into HI (/(r, u) £ H for any u £ H) 

(ii) There is a positive constant k such that 

|/(t,0)| < k for all r £ R 

|/(r,ri) - /(r,r 2 )| < k |n - r 2 | for all r X! r 2 G M w and t e R. 

We deduce from [part (ii) of] (A4) the existence of a positive constant c such that 
\f(r,r)\ < c(l + |r|) for any (r,r) e R x R N . The scaled functions p e £ L°°(Q), 
a e = K)i< tJ <A e L°°(Q T ) NxN , r{-,r) £ C(0,T) (for any r £ R N ) and ¥ £ 
L°°(Qt) are defined as follows: 

p £ (x) = p(^j , al^x.t) = aij (X^) i F(-,r)(t) =1 C^> r J and 
6 £ (x,i) = fc|^Jfor(i,()eQ T . 
Let e > be a small parameter. We will make use of the following notation 



(u,v) s — / p e (x)u(x) ■ v(x)dx (u,v £M), 

JQ 
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which defines an inner product in H making it a Hilbert space. The associated 
norm is denoted in the sequel by | • | e and is obviously equivalent to the natural 
norm of H denoted by • | . 

Given uq G H, we are interested in the asymptotic behavior of the sequence of 
velocity field (u £ ) £> q of the following generalized Ladyzhenskaya equation 
(2.1) ' 

fOg. - div (a s Vu e + b £ |V« £ | P ~ 2 V« £ ) + (t*e ■ VK + Vg £ = p e f e (; u £ ) in Q T 

div u e = in Qt 

u e = Q on dQ x (0, T) 

u E (x, 0) = u°(x) in Q. 

We introduce the functionals (for fixed e > 0) 

a e j (t;u,v)= / (a E Vti) ■ Vvdx + / b e |V«| P ~ 2 Vm • Vvdx 

JQ JQ 

and 




defined for u,v,w G W ' P (Q) N . Then the following estimates hold: 

(2.2) |of(f;u,t>)| < HalU \\Vu\\ L2{Q) \\Vv\\ L2[Q) + u 2 HVull^ \\Vv\\ LP{Q) , 

(2.3) a e j(t;v,v) > v \\Vvf LHQ) + Vl ||Vu|| p p(Q) 

for all u, v E Wq ,p (Q) n . Since p > 1 + ^ (hence p > 2) we have L*>(Q) c — > i 2 (Q) 
so that, by the estimate (|2.2[) we infer the existence of an operator A £ (t) : V — >• V' 
such that 

af (i; u, v) = {A £ (t)u, v) for all u, v G V. 

Hence the existence of a bounded operator A £ : L p (0,T;V) L p '(0,T;V) such 
that 

(.4 £ «,u) = / (A e (t)u(t),v(t))dt for all u,w G L p (0,T;V). 
Jo 

Now, dealing with the trilinear functional 6/, we have that [2"2l |3"8] 

bi(u, v,v) = for all u G V and v G W 1,P (Q) JV 
bi{u,u,v) = -bi(u,v,u) for all w GV and v G Wo'^Q)^. 

The following Sobolev embedding holds true [TJ: VF 1 ' P (Q) L r (Q) for ± = ± - i 
if ~ — jj > 0, and for any r > 1 if | — < 0. In view of the choice of p and N, we 
may always find r > 1 such that f + | = 1, and hence, using Holder's inequality, 
we get 

\bi(u,u,v)\ = \bi(u,v,u)\ <c\\u\\ 2 Lr{Q) ||Vv|| iP(Q) (u,ve V). 
We infer from the above inequality the existence of an element B(u) G V' such that 
(B(u), v) = bj(u, u, v) for all u, v 6 V. 
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This defines a bounded operator B : V — > V' verifying the further property that, if 
11 G L'P(Q,T;Y) then B{u) e L p ' (0,T;V). Indeed taking u € LP(0,T;V) we have 
by the Holder's inequality 

(2-4) II^HII^'co^vo ^ IK*)H^'(Q)^ ■ 

But since W l p (Q) '—t L r (Q), there is a positive constant Co independent of m such 
that 




II B ( W )IIlp'(0,T;V) < C 

where we recall that ||-|| stands for the norm in V. We distinguish two situations: 
the case when p > 3 and the case when 1 + < p < 3. If p > 3, then 2p' < p, 
hence the Holder's inequality again gives (noting that #7 > 1) 



2p> 

\ VP' / T \ VP 

H*)ii 2p 'd* < C i / n«(t)ir* . 



so that 

( 2 - 5 ) II B ( m )IIlp'(0,T;V) < c 2 IMlL>(0,T;V) • 

Now, if 1 + j|^2 < P < 3, then it is easy to get an inequality similar to (|2.5I) for 
N = 2. So we shall only deal with the case N = 3. With this in mind, returning 
to the inequality (|2.4[) which is true for r — 2p' , we have 

/ T f \2/p 

ii b niIl,'(o,t;vo < c 2 1 y n ii«(t)iifv ( g) # J • 

But [2D] (see also [H Lemma 4.1]) 

(2.6) |M| iV(Q ) < C \\u\\ a wl . P{Q) \\u\fe" Q) for any u G I^(Q) 

where a = 3/(5p - 6). But if u G L p (0,T;V) n L°°(Q, T; H), the function t i-> 
||u(i)|| 2a |«(i)| 2 ~ 2a lies in L r °(0,T) for r > 1 satisfying i = ^ + i, and more 
generally, in any Z7(0,T) for 1 < r < Since p' < ^ (for 3 > p > ^-) we 
have that it belongs to L p (0, T). Using Holder's inequality in (|2.6p (with exponent 
r = tt^-t) we have 

2ap' J_ 

II^WII^(o,T ;V o^ C (/ T |l U ^ll Pd *) P Wt)\ i2 ~ 2a)P ' r ' dtj 

(where r' = r/(r — 1)), or equivalently, 

(2-7) II^(' U )IIlp'(0,T;V') — ^ H M Hlp(0,T;V) ll U lli°°(0,T;H) 

for any it 6 L p (0,T; V) D L°°(0,T;H). 

We have proven the following inequalities: 

Ifp>3, ||B(«)|| L ^ (0iT .v') < C2 II«IIL(o,t ; v) forall«G^(0,T;V); 
(2.8) ^ If 1 + ^ <P<3, ll-B^ll^^^.y^ < C||m|| 2 p (0T;V) ||u|fc 2 ( o, T; H) 
for all m G L p (0,T;V) R L°°(0, T; H). 
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The above inequalities will be very useful in the sequel. For the sake of com- 
pleteness, we choose it G H. We are therefore concerned with the existence of a 
solution of (|2.ip . The first result of the work is the following 

Theorem 1. Let 1 + < p < oo. Suppose u° G H. Under assumptions 

(A1)-(A4), there exists (for each fixed e > 0) a couple (u £ ,q £ ) G i p (0,T;V)n 
L°°(0,T;H) x W-^iO^iLP'iQ)) solution to fl2H]). The function u e also belongs 
to C([0, T]; H) and q £ is unique up to a constant function of x: Jq q £ dx = 0. 

Proof. Multiplying Eq. (|2.1[) by v G V and integrating over (0, t) x Q we get 

{uS),v) e + J* {A £ (s)u £ (s) + B(u £ (s)),v) ds = (u°, v) s 
+ f Q (/ (^-, u £ (s)) ,v) e ds for all v G V and a.e. < t < T, 

which, in view of the properties of the operators A £ and B (see especially (|2.8[l ). 
amounts to find a function u £ G L p (0,T;V) such that 

(ill I 

(2.9) —l+A e u E + B(u £ ) = f £ (-,u £ ) in L p (0,T;V), « e (0) = u° in H. 

dt 

Conversely, a solution of (|2.9p will satisfy (|2.ip for a suitable choice of (? e which shall 
be specified later. Therefore, in view of the properties of the operators A e and B 
(see once again (|2.8l0 . we can argue as in [51 Theorem 2.1] (see also [5]) to get the 
existence of a solution to (HH]) in the space L p (0,T; V) CI L°°(0,T;V). It is to be 
noted that the solution u £ satisfies '^f- G LP (0,T;V), so that, by a well-known 
result, u e eC([0,T];H). 

For the existence of the pressure, we have p e f £ (-,u £ ) G L 2 (Qt) n , so that the 
necessary condition of [3SJ Section 4] for the existence of the pressure is satisfied. 
Now, coming back to (|2.ip and denoting there 

w £ = p £ f £ (-, u £ ) -p £ ^+ div {a £ Vu £ + b £ |Vw £ | p ~ 2 Vm e ) - (u e ■ V)u £ 

we have (w £ ,v) = for all v G V (that is for all v G C£°(Q) N with divv = 0) 
where (•, •} is the duality pairing between V(Q) N and T>(Q) N . Next, arguing as in 
the proof of [33 Proposition 5] we are led to w £ G VF _1 '°°(0, T; W -1 ' 13 ' (Q) N ), so 
that there exists a unique g e G W~ 1,oo (0, T; L p (Q)) such that 

Vq £ = w £ , / <7 e da; = 0. 

This completes the proof. □ 
2.2. A priori estimates and compactness. The following result holds. 

Lemma 1. Let 1 + ^ 2 < p < oo. Under assumptions (Al)-(A^) we have the 
following estimates: 

(2.10) sup \u £ (t)\ 2 < C, 

0<t<T 

(2.11) f T \\u £ (t)\\Pdt < C 

Jo 

where C is a positive constant which does not depends on e. 
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Proof. The variational formulation of ([2.1J gives, for any tieV, 

(2.12) j t (u £ (t),v) E + a e I {t-u £ {t),v) + b I {u £ (t),u £ (t),v) = (j (j^,u s (t)j ,v 

In all what follows, C is a generic constant that may change from line to line. 

Taking the particular v — u £ (t) in (|2.12|) and using the relation bi(u £ (t), u £ (t), u £ (t)) = 
0, we get (after integrating over [0,t]) for all t £ [0, T], 

\u e (t)\l + 2v / \\u £ (s)\\ 2 H i, Q)N ds + 2vi \\u £ (s)\\ p ds 
Jo Jo 

< K|J + 2 / /(±, Us ( s )) Ma^ds. 

JO V£ 7 E 

From assumption (A4) we have that 

f(^,u £ (s)) <C(l + |u e ( fl )| B ), 

2/ / (4, «=(*)) l«e(*)| e ds < ^/(l + M*)l e )Ma)l e *» 
Jo V£ ' e Jo 



hence 



< C + C 



[ \u e (s)\ 2 e ds. 
Jo 



We therefore get 

<c + c7 V( s )|^ s 

since sup e>0 I it 1 2 < 00. We readily deduce from (|2.13l) 



M*)l*<C + <? [ \u e (s)\ 2 e ds. 
Jo 



By the application of Gronwall's inequality we find that 
(2.14) \u s (t)\ 2 e < C for < t < T, e > 

where C is a positive constant independent of e and i. But |it e (i)|^ > A -1 |it e (£) | 2 , 
so that (f2~TU|) comes from (j2~T4l) . We also infer from (j2~TB"j) that (f2~TT|) holds true. 
Still from (|2.13p it emerges that 

(2.15) 



sup I ||u e (t)||^L (Q)J v dt < C. 
The next result will be of great interest in the homogenization process. 



□ 



Proposition 1. The sequence (u e ) e >o is relatively compact in the space L 2 (Qx) N ■ 
Proof. Let 

u 6 L*(0, T; V) : — G ^ (0, T; V + L 2 (Q) N ] 



is compact and H ^ L (Q) , we have that 



L 2 (Q) N is corn- 



Since 

pact. We also have L 2 (Q) N » V + L 2 (Q) N , hence the compactness of W in 
L 2 (0,T; L 2 (Q) N ) = L 2 (Q T ) N by a well-known result; see e.g. [22J p. 58, Theorem 
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5.1]. With this in mind we need to check that the sequence (it e ) e >o is bounded 
in W. But since V ^ V + L 2 {Q) N we just need to verify (after (|2XL|) ) that 
(du £ /dt) e> o is bounded in LP (0,T;V). But this comes from the combination of 
(j2lu| -(|CT ]) with (HU and (j235)l . □ 

We can now deal with the estimation of the pressure. Before we can do that, let 
us recall the definition and properties of the Bogovskii operator. Let 



L P (Q) = 6 L P (Q) : J vdx = j 



for 1 < p < oo. We have the following result. 

Lemma 2 ([571 Lemma 3.17, p. 169]). Let 1 < p < oo. There exists a linear 
operator B : L P (Q) — > W ' P (Q) N with the following properties: 

(i) div#(/) = / a.e. in Q for any f £ L P (Q) 

(ii) \\B(f)\\ w u P{Q)N <c( P ,Q)\\f\\ LP{Q) 

(iii) Lf f — div g with g £ L r (Q) N and g ■ n = on <9Q /or some 1 < r < oo 
where n is an outward unit vector normal to dQ, then 

\\B(f)\\ Lr{Q)N <c(r,Q)\\g\\ Lr{Q)N . 

(iv) // / £ C °°(Q) n L p Q (Q) then B(f) £ C °° (Q) w 
With this in mind, the following result holds. 

Lemma 3. Let l+-^pj < p < oo andp' its conjugate. We have q e £ L p (0, T; Lq (Q)) 
with 

(2-16) sup\\q £ \\ LP , (QT) <C 

e>0 

for some positive constant C independent of s. 

Proof. First we know that q £ £ V'(Q T ). Let g £ C^°(Qt) H L p (0, T; Lq(Q)) (recall 
that the dual of Lq (Q) is L P (Q)). In view of [part (iv) of] Lemma [5] let v £ 
C™{Qt) N be such that divv = 5. We have 

(Vq e ,v) = - {q £ ,divv) = - (q £ ,g) , 

that is (because of (|2.1|Q 

(<7e,.g) = - (Vg E ,u) = (-w e ,v) 

where 

u> £ = p e / £ (-, M e ) - + div (« £ Vii £ + b £ \Vu £ \ p ~ 2 V« £ ) - {u e • V)tt £ . 

But we infer from Lemma [T] that 

|((fe,5>| < C (\\v\\ LPiQtT . w i,P (Q} N } + \\v\\ lp{Qt) nJ 

< C ll t 'llLP(0,T;H' o 1 ' ! '(Q)«) 

< C||5IIlp(q t ) > 

the last inequality above being due to [part (ii) of] Lemma We therefore deduce 
from the above inequality that q £ £ L p (Qt), that is q £ £ L p (0, T; Lq (Q)) with 

II&IIl-' Wt) <c 

for a constant C > independent of e. □ 
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3. The ^-convergence method 

We begin this section by collecting some useful tools about almost periodicity. 

3.1. Almost periodic functions. The concept of almost periodic functions is 
well known in the literature. Following [291 Section 2], we present in this section 
some basic facts about it, which will be used throughout the paper. For a general 
presentation and an efficient treatment of this concept, we refer to [5], [3] and |21) . 

Let B(R N ) denote the Banach algebra of bounded continuous (complex- valued) 
functions on R N endowed with the sup norm topology. 

A mapping u : R w — >• C is called an almost periodic function, or a Bohr almost 
periodic function, if u G 6(R ) and further the set of all its translates {u(- + a) : 
a G R N } has a compact closure in B(R N ). We denote by AP(R W ) the set of all 
continuous almost periodic functions on R . AP(1& N ) is a commutative C*-algebra 
with identity. Next, let us denote by Trig(R Ar ) the algebra of all trigonometric 
polynomials, i.e. all finite sums of the form 

u (y) = ak ex p(^fc • y)' £fc e rJV ' a k e c - 

Then ^(t™) C AP(R N ) and further AP(R N ) coincides with the closure of 
Trig(M JV ) in B(R N ). From the above definition, one easily sees that every ele- 
ment of AP(R N ) is uniformly continuous. Moreover it is classically known that 
AP(R N ) enjoys the following properties: 

(P)i u{- + a) G AP(R N ) whenever u G AP(R N ) and for every a G R N ; 
(P) 2 For each u £ AP(R N ) the closed convex hull of {u(-+a)} a&l N in B(R N ) con- 
tains a unique complex constant 371 (it) called the mean value of u, and which 
satisfies the property that the sequence (u e ) e> o (where u E (x) = u{x/e), 
x E R N ) weakly ^-converges in L°°(R N ) to 37t(u) as e -> 0. 

It follows from the above properties that AP(R N ) is an algebra with mean value 
on R N [T^] • Its spectrum is the Bohr compactification of R N , sometimes denoted 
by bR N in the literature and, in order to simplify the notation, we denote it here 
just by /C. The set K. is a compact topological Abelian group. The Haar measure 
on /C is denoted by (3. In view of the Gelfand representation theory of C*-algebras 
we have the next result. 

Theorem 2. There exists an isometric ^-isomorphism Q of AP{R N ) onto C(JC) 
such that every element of AP{R N ) is viewed as a restriction to R N of a unique 
element in C(/C). Moreover the mean value 371 defined on AP(R N ) has an integral 
representation in terms of the Haar measure f3 as follows: 

Wl(u)= [ G(u)d(3 for allue AP(R N ). 

JK 

The isometric ^-isomorphism Q of the above theorem is referred to as the Gelfand 
transformation. The image Q{u) of u will very often be denoted by u. 

We introduce the space AP°°(R N ) = {u G AP(R W ) : D*u G AP(R W ) for every 

a = (ai, . . . , ajv) G N N } where Dy = ay *N ■ F° r m e N (the non negative 

integers) and for u E AP°°(R N ) we set |||u||| m = sup| Q | <m sup ygR ]v \Dyu\ (which 
defines a norm on AP°°(R N )). Then AP°°(R N ) is a Frechet space with respect to 
the natural topology of projective limit, defined by the increasing family of norms 
IH-IIL(mGN). 
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Next, let B AP (R N ) (1 < p < oo) denote the space of Besicovitch almost periodic 
functions on M. N , that is the closure of AP(R N ) with respect to the Besicovitch 
scminorm 

IMI = (limsup— / \u(y)fdy 

V r^+oo \D r \ JB r 

where B r is the open ball of 1* of radius r centered at the origin. It is known that 
B Ap (R N ) is a complete seminormed vector space verifying B q Ap {U. N ) C B p Ap {R N ) 
for 1 < p < q < oo. Using this last property one may naturally define the space 
B AP (R N ) as follows: 

Bf P (R N ) = {/ G n 1 < p<oc i?^ P (R A ') : sup ||/|| p <oo}. 

l<p<oo 

We endow B Ap (R N ) with the seminorm [f]^ = sup 1<p<oc ||/|| p , which makes it 
a complete seminormed space. We recall that the spaces ^^(R^) (1 < p < oo) 
are not Frechet spaces since they are not separated. The following properties are 
worth noticing [2%l [28] : 

(1) The Gelfand transformation Q : AP(R N ) — >• C(JC) extends by continuity to 
a unique continuous linear mapping, still denoted by Q, of B Ap (R N ) into 
L P (K,), which in turn induces an isometric isomorphism Qi, of B AP (R N )/M = 
B P AP (R N ) onto LP(IC) (where N = {u S B p Ap (R N ) : Q(u) = 0}). More- 
over if u G B AP (R N ) n L°°(R N ) then Q{u) G L°°{K) and \\Q{u)\\ LX{K) < 

\\ u \\l°°(r n )- 

(2) The mean value 9Jt, defined on AP(R N ), extends by continuity to a positive 
continuous linear form (still denoted by SUt) on ^^(R^) satisfying ?0t(u) = 
J K G{u)d(3 and Wl(u(- + a)) = M(u) for each u G B AP (R N ) and all a G R N , 
where u(- + a)(y) = u(y + a) for almost all y G R N . Moreover for u G 
B AP (R N ) we have ||u|| p = [M(\u\ p )] 1/p . 

Spaces of almost periodic functions with values in a Banach space are defined 
in a natural way, we refer to [4] for details. Keep the following notations in mind: 
AP(R N ; C) = AP(R N ) and B p Ap (R N ;C) = B p Ap {R N ). 

We can also define the notion of almost periodic distributions. To do this, let 
B oc (R N ) denote the space of all C°° functions in R N that are bounded together 
with all their derivatives of any order. We equip Boo(R N ) with the locally convex 
topology defined by the increasingly filtered separating family of norms || |-| \\ m (m £ 
N). As usual, let ^(R^) denote the subspace of ^^(R^) consisting of functions 
with compact support. We denote by ^^(R^) the topological dual of B 00 (R N ). Wc 
recall that B'^R™) is not a subspace of V(R N ) (the usual space of distributions 
on R N ) since V(R N ) is not dense in S OQ (M JV ); see [33]. We endow ^(K^) with 
the strong dual topology. For a distribution T G V(R N ) we define its translate 
r a T {aeR N ) as follows: 

(r a T, <p) = (T, tp(. - a)) for any <p G V{R N ). 

With this in mind, we say that a distribution T is an almost periodic distribution 
if T e B' oa (R N ) and further the set of translates {r Q T : a E R N } is relatively 
compact in B' oa (R N ). We denote by B' AP (R N ) the space of all almost periodic 
distributions on R N . B' AP (R N ) is endowed with the relative topology on 23^ (R w ). 
The following result holds [33] p. 206, Section 9]. 



/ / 
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Proposition 2. The following assertions are equivalent. 

(i) T G B' AP (R N ); 

(ii) T G B^R^) is a finite sum of derivatives of functions in AP(R N ); 

(iii) T G B^(M W ) and T * tp E AP°°(R N ) for any ip £ V(R N ). 

It can be easily checked that <^T G B^pQR^) whenever tp G AP°°{R N ) and 
T G B' AP (R N ). It follows from part (iii) of the above proposition that AP°°(R N ) is 
dense in B' AP (R N ). Since B\ P (R N ) is the closure of AP°°(SL N ) in the Besicovitch 
seminorm whereas B' AP (R N ) is the closure of AP°°(R N ) in B' 00 (R N ) it follows that 
B AP (R N ) is continuously embedded in B' AP (R N ). As a first consequence of this, 
we have the 

Lemma 4. The mean value 9Jt considered as defined on AP°°(R N ) extends by 
continuity to a unique continuous linear form still denoted by 971, on B AP (R N ). 

Proof. The mapping 9JT being continuous on AP°°(R N ) (endowed this time with 
the relative topology on B^R^)), the result is therefore an immediate consequence 
of the density of AP°° (R N ) in B' AP (R N ). □ 

OJt(T) is still called the mean value of T G B' AP (R N ). Now, let T G B^ P (R Ar ). 
Using the preceding lemma and the density of AP°°(R N ) in B' AP (R N ), we may 
define the spectrum of T as Sp(T) = {k G R w : 9#(7 fe P) ^ 0}, a countable subset 
of R*, where 7j.(y) = exp(ifc • y) for y G R w and 7 fc stands for the complex 
conjugation of j k . We know from [33l p.208] that if M{^j k T) = for any k eSp(T) 
then T = 0. 

Let R^+ x = x 1 T denote the space I ff xl with generic variables (y,r). 
It is known that AP(R^+ X ) = AP(R T ; AP(R^)) is the closure in B(R^+ 1 ) of the 
tensor product AP(R^) ® AP(R r ) [II]. We may sometimes set A y = AP(R^), 
A T = AP(R T ) and A = AP(R^+ X ). Correspondingly, wc will denote the mean 
value on A$ (( = y, r) by 9Jtf . 

Now let 1 < p < oo and consider the A-parameter group of isometries {T(y) : 
y G R w } defined by 

Tfo) : ^ P (M W ) -> B^ P (E JV ), T(y)(u+Af) = u(- + y) + Af for u G B^ P (1 W ). 

Since AP(R Ar ) consists of uniformly continuous functions, {T(y) : y G R^} is a 
strongly continuous group in the following sense: T(y)(u+Af) — > u+Af in B AP (R N ) 
as \y\ — > 0. Using the isometric isomorphism Q\ we associated to {T(y) : y G R N } 
the following TV-parameter group {T(y) : y G R w } defined by 

T(y) : LP(X:) -> L"(K) 

T(y)£i (u + AO = Gi(T(y)(u + Af)) = Qi(u(- + y) + AT) for u £ B' Ap (R N ). 

The group {T(y) : y £ R N } is also strongly continuous. The infinitesimal generator 
of T{y) (resp. T(y)) along the ith coordinate direction, denoted by Di tP (resp. di. p ) 
is defined by 

Di iP u = lim s _>o (T(sei)u — u) in B AP (R N ) 
(resp. di tP v — lim s _j.o s _1 \T(sei)v — v) in L' P (1C)), 

where we have used the same letter u to denote the equivalence class of an element 
u G B P AP (R N ) in B P AP {R N ), e, = (%)i<j<at (<% being the Kronecker 6). The 
domain of Di tP (resp. in B AP (R N ) (resp. L V {JC)) is denoted by X>j iP (resp. 
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Wi.p). It emerges from the general theory of semigroups [HI Chap. VIII, Section 
1] that the following result holds. 

Proposition 3. X>, )P (resp. Wi iP ) is a vector subspace of B AP (R N ) (resp. L P (IC)); 
Di tP : T>i >p — > B AP (R N ) (resp. di. p : Wi, p — > L P (K,)) is a linear operator; T>i jP (resp. 
Wi.p) is dense in B P AP (R N ) (resp. L P (JC)). and the graph of Di iP (resp. di tP ) is 
closed in B P AP {R N ) x B P AP (R N ) (resp. W{K) x L P {K)). 

In the sequel we denote by g the canonical mapping of B Ap (R N ) onto B AP (R N ), 
that is, g(u) = u + N for u G B AP (R N ). The following properties are immediate. 
The verification can be found either in [39l Chap. Bl] or in [4]. 

Lemma 5. Let 1 < i < N. (1) If u G AP 1 (R N ) then g{u) G V i<p and 
(3.1) D hp g(u) ^gij^- 



dyi 

(2) If u G V itP then Gi(u) G Wi iP and Qi{D itP u) — d i)P Qi{u). 

We define the higher order derivatives as follows: Dp* = Df p o ■ • • o D^f (resp. 
d p = dl p o • • • o for a = (oi, ...,a N ) G N w with £>£ = D hp o ■ ■'■ o D hp , 

ai-times. Now, set 

B^QR W ) = f)f =1 V hP = {»e B P AP {R N ) : Di, p u G B^R^) VI < » < N} 

and 

Pap(1") = {u G B^CR^) : £>;> G B^p(lR JV ) Va G N^}. 

One can show that V AP (R N ) is dense in B Ap (R N ), 1 < p < oo. Equipped with 
the norm 

/ AT \ VP 

h\\ B teP») = (iK+EiiA^in ( U e^(M w )), 

B^'p (IR^) is a Banach space: this comes from the fact that the graph of Di p is 
closed. 

The counter-part of the above properties also holds with 

W^iK) = nf =1 W i<p in place of B A P P (R N ) 

and 

V{K) = {u G L°°(/C) : d^u G L°°{1C) Va G N N } in that of V AP (R N ). 

Moreover the restriction of Gi to B AP (R N ) is an isometric isomorphism of B AP (R N ) 
onto W 1,p (/C); this is a consequence of [Part (2) of] Lemma[5j 
Let u G T>i^ p (p > 1, 1 < i < N). Then the inequality 



,-i 



s 



(T{sei)u — u) — Di. p u ][ ^ l < c ||s l {T(sei)u — u) — Di p i 



for a positive constant c independent of u and s yields Di^u — Di. p u, so that Di p 
is the restriction to B Ap (R N ) of D^i. Thus for all u G 2\oo we have it G 2\ p 
(p > 1) and Di^u = Di. p u for any 1 < i < N . The next result will be useful in 
the sequel. 

Lemma 6 ([2H1 Lemma 2]). We have V AP (R N ) = g(AP°°(R N )). 



HOMOGENIZATION OF LADYZHENSKAYA EQUATIONS 



13 



From now on, we write u either for Q(u) if u g B Ap (R N ) or for Q\(u) if u g 
B AP (R N ). The following properties are easily verified (see once again either [39l 
Chap. Bl] or 0]). 

Proposition 4. T/ie following assertions hold. 

(i) = for all u g V AP (R N ) and a g N N ; 

(ii) J^. di tP ud[3 = for all u g 2?i !?) and 1 < i < N ; 

(iii) Di tP (u(f>) = uDi^tf) + (f>Di tP u for all ((/>, u) g X'ap(K A ') x 2?i,p flfwi 1 < * < 

n' 

From (iii) in Proposition |4] we have 

/ 0d itP ud0 = - ud i<00 (j)df3 for all (u,<f>) g £>i, p x T> AP (R N ). 
Jk. Jk 

This suggests us to define the notion of distributions on D AP (R N ) and of a weak 
derivative. First, let us endow V AP (R N ) = g(AP°° (R N )) with its natural topology 
defined by the family of norms N n (u) = sup| Q | <n sup ygR w \D^ u(y)\, integers n > 0. 
So topologized, T> AP (R N ) is a Frechet space. We denote by T> AP (R N ) the topolog- 
ical dual of V AP (R N ) and endow it with the strong dual topology. The elements 
of V AP (R N ) are called the distributions on T> AP (R N ). There exists a relationship 
between the space of almost periodic distributions B' AP (R N ) and V' Ap (R N ). To be 
more precise, we have the following result which allows us to view B' AP (U. N ) as a 
proper subspace of V Ap (R N ). 

Proposition 5. Given any T <E B' AP (R N ), the linear form $(T) : V AP (R N ) -)■ C 
given by 

($(T), Qi<p)) = SK^T) g AP»(K W )) 

zs a distribution on T> AP (R N ). The map $ : T i— > $(T) is a linear continuous 
embedding of B' AP {R N ) into V' Ap (R N ). 

Proof. Let (f?(¥>j))j be a sequence of functions in 2? j 4p(R a ') such that f?(y,-) — > 
in V AP (R N ); then ^ -> in AP°° (R N ). In fact assuming gfa) -> in 2? j4 p(M Ar ) 
leads to QiiQiD^ipJ)) = G{D%tp 5 ) -> in C(/C) for any a g N N , so that, as 
||£>™^-|| = VjOHoo (see Theorem[2]) our claim is justified. It follows that 

tfj -> in S^ P (R JV ), hence ->• in B' AP {R N ). We deduce from Lemma H 
that M(tpjT) 0, from which the continuity of $(T) on P AP (R Ar ). Moreover, 
if $(T) = then £E%>T) = for all tp g AP 00 ^), hence SDt(7 fc T) = for any 
k g R N . This yields T = (see [331 P-208]). $ is then injective. It remains to 
check that <& is continuous. But this follows at once by the continuity of the mean 
value Wl on B' AP {R N ) (see Lemma SJ). □ 

Remark 1. The map $ defined above is not surjective since B' AP {R N ) is not the 
topological dual of ^"(R-^) (see [33J p. 207]). Since it is injective and continuous, 
B' AP (R N ) can be viewed as a proper subspace of V AP (R N ). 

The weak derivative of u G V AP (R N ) is defined as follows: for any a g N , Z? Q u 
stands for the distribution defined by the formula 

(L> Q ii,</>) = (-l)H (u,D^(f>) for all g V AP {R N ). 
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The space V AP (R N ) being dense in B AP (R N ) (1 < p < oo), it follows that 
B AP (R N ) C T>' Ap (R N ) with continuous embedding, so that the weak derivative 
of any / G B Ap (R N ) is well defined and verifies the following functional equation: 

(D a f, 4>) = (-1) H / fd^dp for all 4> G V AP (R N ). 
Jk 

As a special case, for / G T> i>p we have 

- / fdijdp = [ WiJW Vcf> G V AP (R N ). 
Jk Jk 

Hence we may identify D ip f with D ai f, a.^ = (5ij)\<j<N- Conversely, if / G 
B AP (R N ) is such that there exists ft G B P AP (R N ) with (_D Qi /, <j>) = - f K Mdfi for 
all (f> G P^pCR^), then / G 2\ p and Di, P / = fi- This allows us to justify the fact 
that B AP (R N ) is a Banach space under the norm |H| e ^p( R iv). The same is true for 

W^iK). Moreover V AP (R N ) (rcsp. £>(£)) is a dense subspace of B AP (R N ) (resp. 

We now define the appropriate space of correctors. For that, let p G B 2 AP (R N ) n 
L°°(R N ) be freely fixed with DJt(p) > (OTt(p) the mean value of p). We begin by 
defining the following space: 

B'/pJ^) = {« e : = 0}. 

We endow it with the seminorm 

/ n \ Vp 

NI#,p= (EH^IIpJ («e^ p (K^)). 

Since 9Jt(p) ^ it follows that |H|# p is actually a norm on B AP JR N ). Under 

this norm B AP p (R N ) is unfortunately not complete. We denote by B^ Ap (R N ) its 
completion with respect to the above norm and by J p the canonical embedding 
of B APp (R N ) into B*f AP (R N ). It can be easily checked that V APp (R N ) = {u G 
V AP (R N ) : Wl(pu) = 0} is dense in B APp (R N ). It holds that: 

(Pi) The gradient operator D p = (D 1>p , D NjP ) : B APp (R N ) -> (B P AP (R N )) N 
extends by continuity to a unique mapping D p : B^ Ap (R N ) — > (B AP (R N )) N 
with the properties 

Di p — Di^p o Jp 

and 

\ i/p 

II"H#,p = (EH^IIp iovueB'/ AP {R N ). 




(P 2 ) The space J p (B AP p {R N )) (and hence J p {V APp (R N ))) is dense in B^ p Ap (R N ). 



Moreover the mapping D p is an isometric embedding of Bj. AP (R N ) onto a closed 
subspace of (B AP (R N )) N , so that B 1 / AP (R N ) is a reflexive Banach space. By 
duality we define the divergence operator dhy : (B P AP (R N )) N — > (£^ P (IR W ))' 
(p'=p/(p-l))by 

(3.2) (div,/u,t;) - - (u,D p v) for w G B^ Ap (R N ) and u = («<) G (B P AP (R N )) N , 
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where (u, D p v} = X)j=i Ik. Uidi tP vdf3 . The just defined div p ' operator extends the 
natural divergence operator defined in T>ap(R- N ) since Di p f = Di tP (J p f) for all 

feV AP (R N ). 

Now taking in (|3.2[) u — D p >w with w 6 B Ap (R N ) being such that D p 'W G 

(23^(1^))^, we define the Laplacian operator on B AP (R N ) (denoted here by A p >) 
as follows: 

(3.3) (A p ,w,v) = (div p >(D p ,w),v) = - (D p ,w 7 D p v) for all v G B^ P Ap (R N ). 

If in addition v = J p (4>) with 4> G V A p(R N )/C then (A p <w, J P {4>)) = - (D p >w, D p cf)), 
so that, for p = 2, we get 

(A 2 u>, J 2 (0)) = (w, A 2 <j)) for all w G ^(R^) and <£ G ^Ap(R Ar )/C. 

The following result is also immediate. 

Proposition 6. for u G AP°°(R JV ) we have 

A p g(u) = g(Ayu) 

where A y stands for the usual Laplacian operator on M.^ . 

We end this subsection with some notations. Let / G B AP (R N ). We know that 
D ai f exists (in the sense of distributions) and that D ai f = Di. p f if / G So 
we can drop the subscript p and henceforth denote Di p (resp. di iP ) by djdiji (resp. 
d{). Thus, D y will stand for the gradient operator (d/dyi)i<i<N and divj, for the 
divergence operator div p . We will also denote the operator D i-p by d/dyi. Since 
J p is an embedding, this allows us to view B A ^ p (R N ) (and hence V AP>p (R N )) as a 
dense subspace of Bj AP (M. N ). Di iP will therefore be seen as the restriction of Di. p 
to B2 Pjp {R N ). Thus we will henceforth omit J p in the notation if it is understood 
from the context and there is no risk of confusion. This will lead to the notation 
D p = D y = (d/dy.i)i<i<N and d p = d — (di)i<i<N- Finally, we will denote the 
Laplacian operator on B AP (R N ) by A y . 

3.2. The ^-convergence. Let A y = AP(Ry ) and A T = AP(R T ). We know that 
A = AP(M^+ 1 ) is the closure in B(R y ^+ 1 ) of the tensor product A y ® A T . We 
denote by K, v (resp. /C T , K.) the spectrum of A y (resp. A Tl A). The same letter 
Q will denote the Gelfand transformation on A y , A T and A, as well. Points in ]C y 
(resp. K. T ) are denoted by s (resp. sq). The Haar measure on the compact group 
JCy (resp. fC T ) is denoted by fi y (resp. j3 T ). We have /C = JC y x fC T (Cartesian 
product) and the Haar measure on K, is precisely the product measure (3 = P y ® /3 T ; 
the last equality follows in an obvious way by the density of A y ® A T in A and by 
the Fubini's theorem. Finally, the letter E will throughout denote (£ n )neN with 
< e n < 1 and such that e„ — > as n — > oo. In what follows, we use the same 
notation as in the preceding section. 

Definition 1. A sequence (it e ) e >o C L p (Qx) (1 < p < oo) is said to weakly X- 
converge in L P (Q T ) to some u G L P (Q T ; B AP (M.y J: 1 )) if as e — » 0, we have 

J QT u £ (x,t)f(x,t,^,^)dxdt 

-> 5$Q TyK u (x,t,s,S())f(x,t,s,s )dxdtdl3 

for every / G L p ' (Q T ; A) (1/pf = 1 - 1/p), where u = 5i ou and / = Q x o (gof) = 
Q o f. We express this by writing u e — > uq in L p {QT)-weak E. 
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We recall some important results whose proofs can be founded in [24l Theorems 
3.1 and 3.5] (see also [28l Theorems 3.1 and 3.6]). 

Theorem 3. Let 1 < p < oo. Let {u £ ) £ ^e be a bounded sequence in L p (Qt). Then 
there exists a subsequence E' from E such that the sequence (u £ ) £ ^e' * s weakly 
^-convergent in L p (Qt). 

The next result is of capital interest in the homogenization process. 

Theorem 4. Let 1 < p < oo. Let (u £ ) £e E be a bounded sequence in L p (0, T; W ,P (Q)). 
Then there exist a subsequence E' of E and a couple of functions 

K, ui) G £ p (0, T; W^ P (Q)) x U>(Q T ; B P AP (R T ; S^ P «))) 

such that, as E' 3 e — > 0, 

(3.5) u £ ->■ u in L p (0, T; W^ p (Q))-weak; 

(3.6) ^ -> ^ + ^ m L p {Q T )-weak E, 1 < * < N. 
dx t axi dyi 

We shall also deal with the product of sequences. In this respect, we give a 
further 

Definition 2. A sequence {u £ ) £>0 C L P (Q T ) (1 < p < oo) is said to strongly E- 
converge in L P {Q T ) to some u G L P (Q T ; B Ap (RyJ: 1 )) if it is weakly ^-convergent 
towards uq and further satisfies the following condition: 

( 3 - 7 ) \K\\ L p(Q t ) -> INI| L p(Q t x/C) • 

We denote this by writing u £ — > uo in L p (QT)-strong E. 

This being so we have the following. 

Theorem 5 ( [31 , Theorem 6]). Let 1 < p, q < oo and r > 1 be such that 1/r — 
l/p+ 1/q < 1. Assume (u £ ) £ ^e C L q (Qx) is weakly "E-convergent in L q {Qx) to 
some u 6 L q (QT] B q Ap {M.y , and (v £ ) £e E C L P (Q P ) is strongly ^-convergent in 
L p (Qt) to some Vq 6 L p (Qt; B AP (M.y^ 1 )). Then the sequence (u £ v £ ) ££ e is weakly 
Yi-convergent in L t (Qt) to uqVq. 

As a consequence of the above theorem the following holds. 

Corollary 1. Let (u £ ) £eE C L P {Q T ) and (v £ ) eeE C L p ' (Q T ) H L°°(Q T ) (K P < 
oo and p' = p/(p — 1)) be two sequences such that: 

(i) u £ — > uo to L p {QT)-weak E; 

(ii) u e — >■ wo in- i p (Qr)-strong E; 

(iii) (w e ) ee £ is bounded in L°°(Qt)- 

Then u £ v £ — > itoVo in L p (Qx)-weak E. 

Regarding the reiterated E-convergence we refer the reader to [T3J ST] . 
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4. HOMOGENIZATION RESULTS 

Throughout this section we make the following assumption on the functions 
involved in (|2.1[) : 

(A5) Almost periodicity. We assume that the functions Oy and b lie in 

B% P (R^+ 1 ) n L°°(M^+ 1 ) for all 1 < i,j < N. We also assume that the 
function p belongs to B AP (R N ) n L°°(R N ) with M y (p) > 0. Finally the 
function r H> /(r, r) lies in (AP(R T )) N for any r € R N . 

4.1. Preliminary results. Let 1 + < p < oo. It is a fact that the topological 
dual of B AP (R T ;B^ P AP (M^)) is B% P (R T ; [B^ P AP (R*)]'); this can be easily seen from 
the fact that B^ P AP (R*) is reflexive (see Section 2) and B AP (R T ] B^ P AP (R*)) is 
isometrically isomorphic to L P (IC T ; B^ Ap (Ry)). We denote by (, } (resp. [, ]) the 
duality pairing between S^p(R^) (resp. B AP (R? ; B# AP (Ry ) ) ) and [B^ P AP (R^)}' 

(resp. B^( R t! [B^ P (lJ')]')). Hence, for u e B P AP (R T ; [B^ P Ap(K )]') and w e 
S* p (R T ;B^ p (Rf)), 

[u,v]= {u(s ),v(s ))dP T {s ). 

For a function ^> € P j 4p iP (]R w ) we know that -0 expresses as follows: ^ = Qyirfi) 
with 0! 6 AP p °°(R^) = {v € AP°°(R^) : SUl(pVi) = 0}, where g y denotes the 
canonical mapping of B Ap (Ry) onto £5^(11^); see Section 2. We will refer to tp 1 
as the representative of ip in AP 00 (R^ ) . Likewise we define the representative of 
i> e V AP {R T ) <g) V A p,p(R N ) as an element of AP°°(R T ) ® AP p °°(R^) satisfying a 
similar property. 

With all this in mind, we have the following 

Lemma 7. Let tp e C^°(Qt) <8> [2?ap(Rt) ® 2?AP,p(R^)] a»i<f Vh ^ e representa- 
tive in C^{Q T ) ® [AP°°(R T ) ® AP p °°(R^)]. Let (u e ) £e£ ;, P' and (u ,ui) &e as in 
Theorem [4] Then, as E' 3 e — >• 0, 

—u e p e ip\dxdt^ I {pui(x,t) 7 ip(x 7 t)} dxdt. 



Qt t JQt 
Proof. We recall that for V>i as above, we have 

ip\{x,t) = ^M, ~, ^ for (x,t) S Q T - 

This being so, since pViOM, G ^Pp(R^) = G AP(R^) : £QTj,(/ttt) = 0}, 
there exists Q]5] a unique </> G C£°(Qt) <8> [AP°°(R T ) ® AP(R^)] such that pip 1 = 
A y (f>. The result follows at once by the application of [26l Lemma 3.4] (see also 
[T2l IT4"] for some periodic versions of this lemma and their proofs) . □ 

For u £ £>^ P (R r ) we denote by 9/<9t the temporal derivative defined exactly as 
its spatial counterpart djdyi. We also put do — Giid/dr). d/dr and do enjoy the 
same properties as d/dyt (see Section 2). In particular, they are skew adjoint. Now, 
let us view pd/dr as an unbounded operator defined from U = £?^ P (R T ; £>_^ P (R^)) 

into W = B AP (R T ; [B^ p A p( R y)Y)- Similar to Ql pp. 1243-1244], it gives rise to an 
unbounded operator still denoted by pd/dr with the following properties: 
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(P)i The domain of pd/dr is W = {v e U : pdv/dr e U'}; 
(P) 2 pd/dr is skew adjoint, that is, for all u,v € W, 



dv 




du 









(P) 3 The space £ = £>ap(Rt) ® X>Ap P (Rf ) is dense in W. 
The above operator will be useful in the homogenization process. This being so, 
the preceding lemma has an important corollary. 

Corollary 2 (37, Corollary 1]). Let the hypotheses be those of Lemma^\ Assume 
moreover that u% G W. Then, as E' 3 e — > 0, 



di/j e 
eu e p e ——dxdt 
at 



Qi 



p-^-{x,t),ip{x,t) 
or 



dxdt. 



Now assume that ui = (u\) e {B^ AP ) N and t/> = (?/>*) £ (X' AP (R T )(g)DAP, /9 (M^ )) JV - 

Then we may still define [p^r(x, t), tp(x, t)], but this time as J2iLi[p~^r( x ' ^(^j *)]• 
We will use that notation in the sequel. 

We end this subsection with one further result. 



»JV 



Lemma 8 ([371 Lemma 5]). Let h : R 
verifying the following conditions: 

(i) \h(r, n) — h(r, ra)| < k \r\ — r 2 | for any r e 

(ii) h(-,r) e AP(R) /or a/Z r £ R N . 
Let (u e ) e be a sequence in L 2 (Qt) n such that u e 
where Uq € L 2 (Qt) n ■ Then, setting h £ (u e )(x,t) = 
e -> 0, 

h e (u e ) -> h(-,u a ) in L 2 (Q T ) N -weak E 



N be a continuous function 
. and all r\, r% G R N ; 



► ito m L {Qt) as e — > 0, 
h(t/e 2 ,u £ (x,t)) we have, as 



4.2. Homogenization results. Before we can state the homogenization result for 
(|2.1j) we need a few notations. We begin by noting that the space C^°(0, T) eg) V is 
dense in L p (0,T; V). Next, we introduce the space 

B^ = {«e(^ P (tT: dF y u = 0} 

where divyit = YliLi 9u l /dyi, and its smooth counterpart 
Vdiv = {ue (V APtP (R N )) N : divyti = 0} 

= {ue (V AP<p (R N )) N : M(pu) = and div y u = 0}. 
The following result holds. 
Lemma 9. The space Vdw is dense in B^?. 

Proof. This follows exactly in a same way as the proof of [415] Lemma 2.3]. □ 
Now, let 



u e 



^ P (R T ;^):p^£^ P (R T ;(^)')}; 



We set 
and 



-r-oo 
•M) 



Wl' p = LP(0,T;V) x LP(Q T ;W div ) 
[C o °°(0,T) ® V] x [C °°(Q T ) ® (X>ap(Kt) ® Vdiv)] ■ 
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Thanks to Lemma M and to the density of £ = V A p(R T ) ® V A p P (Ry) in W, we 
have the density of J-£° in F ' p . 

Let (u e ) e( zE be a sequence of solution to (I2.1[) (or to (|2.9p Where we assume that 
_B is an ordinary sequence tending to zero with e. In view of Proposition [I] there 
are a subsequence E' of E and a function Uq E L 2 (Qt) n such that, as E' 3 e — > 

(4.1) u E -> w in L\Q T ) N . 

In view of (|2.11[) and by the diagonal process, one can find a subsequence of (u E ) ee E> 
(not relabeled) which weakly converges in L P (Q, T; V) to the function «o (this means 
that ito € L p (0,T;V)). From Theorem [4] we infer the existence of a function 

Api^r] B#ap) 



Ui = (u\)i<k<N G L p {Qt\ £>^ p (R t ; Si^p)^) such that the convergence result 



du F dun dui „ , „ „ v 
(4.2) ^ - ^ + ^ in ^(Q^-weak E (1 < i < N) 

holds when i?' 3 e — > 0. Owing to Lemma [3] (see (|2 . 1 6[) therein) there exist a 

subsequence of E' (still denoted by E') and a function q E L p ' (Q T ; B AP (R N+1 )) 
with J qdx = such that 



Wr recall that ^ = I §^ 



(4.3) g e -> g in L p (Q T )-weak E as E' 3 e -> 0. 

(u„ = (ug) 1<fc<JV ) and &f = (M) 

l<k<N ~ ~ ay * \ ay * J l<k<N 

Now, let us consider the following functionals: 

aj(u,v) = ^2 a i3 (s,s )B j u k B i v k dxdtdl3 



QtXIC 

b(s, s ) \Bu\ p ~ z BiVpMiui dxdtdfi 



N 
i,3=l 



where U>jU k = ^p- + djU k {djU k = Q\ (Jg^j an d the same definition for J}iV k ) and 
Bit = (Djit)i<j<Ar with H)jU = (H>jU k )i<k<N', 



N 

bi(u,v,w) — u -—9-Wodxdt 



JJq t dx, 



,fe=i 



for u = (uq,Ui),v = (vq,Vi),w — (wq,Wi) E ¥q' p . The functionals ai and bi are 
well-defined. Next, associated to these functionals is the variational problem 
(4.4) 

u = (it , wi) e ^o P '■ 

—$Jt y (p) Jq t uo ■ tp' dxdt — Jq t pui,^^ dxdt + aj(it, #) + bj{u,u,&) 
= J Qt m(pf(;u )) ■ xp dxdt for all $ = (Vo.V'i) £ ^o°°- 
The following global homogenization result holds. 
Theorem 6. The couple (iio,i*i) determined by (|4.ip - (|4.2j) solves problem (|4.4|l . 

Proof. From the equality div u £ = we easily derive div ito — and div^iti = 0. 
In order to show that it = (ito,iti) E Fq P we need to verify that du\{x,t)/dr E 
B Ap (M. T ; {B^?)') for a.e. (x,t) E Qt- But this will be done later. Now, we first 
need to show that it solve (|4.4j> . For that, let * = (i/> ,i(>x) with if> E C^(Qt) N 
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and ip x = = {g y {^. j ))i< ] <N e C °°(Q T ) <8> CD AP (R T ) $ Vdiv). Q y being the 

canonical mapping from B p Ap {M.y) into B p AP {M.y); define $ e as follows: 



$ e (x, t) = ip (x, t) + eip Cx, t, j for (x, t) e Q T - 

Then we have <fr 6 G Cq°(Qt) N and, using $ e as a test function in the variational 
formulation of (|2.1[) we get 

rT 

(4.5) 



/ p £ u £ ■ ^-^-dxdt + j a 6 Du £ ■ DQ £ dxdt + / bi(u £l u £ , $ £ )dt 
Jq-t &t JQt Jo 

q £ div & £ dxdt 



b £ \Du, 



IP-2 



Z?u F • DQ F dxdt 



[ ( P s r-(t,u £ (t)),$ £ )dt. 

Jo 



We pass to the limit in (|4.5I) by considering each term separately. First we have 



p u £ ■ — — dxdt — 
dt 

In view of Corollary [5] we have that 
dip 



d ^°-dxdt + e 



dt 



F 7 7 

p u £ ■ - dxdt 



Qt 



dt 



P Us ■ 



Qt 



dt 



-dxdt 



Qi 



pui 



dip, 



dr 



dxdt when E' 3 e -> 0. 



On the other hand, due to (|4.1[) we get 
dip Q 



p u £ ■ 



Qi 



dt 



f dip 
dxdt -> OJty(p) / Mo • —^-dxdt as E 1 3 



Hence, as E' 3 e — > 



Qt 



dt 



dxdt -3- M y (p) [ u ■ ^-dxdt + f 

JQt w JQt 



pui 



dip! 



dr 



dxdt. 



Next, it is an usual well known fact that (see e.g., [T31 Lemma 3.5]), the conver- 
gence result (14.21) together with the weak E-convergence of the sequence (D<& £ ) £ to 
D<&, imply 



N 

a e Du £ ■ D<& £ dxdt -3 

Qt i,j,k=l 



o ii7 -(s, s )D j u k D l v k dxdtd/3. 



Considering the next term, we use the monotonicity property to have 

(4.6) / (b £ \Du £ \ p ~ 2 Du £ — b £ \D& £ \ P ~ 2 D^ £ ) ■ (Du £ — DQ £ )dxdt > 0. 

JQt 

Owing to the estimate (|2.1ip we infer that 

,p-2 



sup 

e>0 



b e \Du F 



Du e 



LP'(Q T ) NXN 



< CO, 



so that, from Theorem[3j there exist a function \ € L p (Qt', B ap (M. n )) NxN and a 
subsequence of E' not relabeled, such that b £ \Du £ \ p ~ 2 Du £ — > x in L p> (Qt) NxN - 
weak £ as E' 3 e — > 0. We therefore argue as in the proof of [551 Theorem 3.1] (see 
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also [HJ Theorem 3.2]) to pass to the limit in (|4.6|) (as E' 3 e -> 0) and get 
(4.7) // (x-b\B$\ p ~ 2 B<I>)) ■ (Bu - B&)dxdtdf3 > 

for any $ G J 7 ^ where, as above, Bit = Z?ito + du\ (u = (uo,Ui)) and B<fr = 
-D'j/'o + S^]^. By the density of J 7 ^ in Fj' p and by a continuity argument, (|4.7[) still 
holds for $ £ F ,p . Hence by taking 3? = u + Xv for i? = (i> , t>i) £ F ' p and A > 
arbitrarily fixed, we get 

(x-b(-,Bu + XTDv)) ■ BvdxdtdPdj3 > for all v G Fg' p . 

QtxK 

Therefore by a mere routine, we deduce that \ = 6 |-D«o + D y u\ \ P 2 (Duq+D v ui). 

The next point to check is to compute the lim e _j.o Jq bi(u e , u £ , <& e )dt. We claim 
that, as e — > 0, 

rT 



/ bi(u £l u £ ,<& £ )dt — ► bi(u,u,&) 
Jo 



Indeed, by the strong convergence result (|4.1[) in conjunction with the Theorem 
our claim is justified. 
We obviously have that 

{ {p £ f{t,U e {t))^ e )dt = f p S f(t,U E (t))-*edxdt 
Jo JQ T 

p e f(t,u E (t)) ■ ip dxdt 



+e p £ f £ (t,u E (t))-ip E dxdt. 
Jq t 

It can easily be seen that the second integral of the right-hand side goes to zero 
with e € E' . For the first one, taking p s ip as a test function and using Lemma [51 
we are led to 

p e / £ (-, u e ) ■ ip dxdt -)• / / p/(-, u ) ■ ip dxdtd(3 

Qt JJQtxIC 



Tt(pf(-,u )) ■ ip dxdt. 
Now, on the basis of (|4.3|) there is no difficulty in showing that 
/ q s div & E dxdt — > // qdiv xp dxdtdp. 

JQt JJQtxIC 

It emerges from the above study that u = (uq,Ui) and go = 90 f t(Q r ) := J^qd/3 
satisfies 

(4 8) "^(^Qr u ° ' ^'odxdt - J Qt pui,^ dxdt +Sj(w,*) + 6, (u, it, *) 
- Jq t qodiv^ dxdt = J Qt Wl(pf(-,u )) ■ ip dxdt 

for all $ = (?/> , ■»/>!) as above. Choosing in particular <I> G we end up with 
(u , Ui) solving (|4.4p since divi/> = in that case. □ 
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Problem (|4.4|) is the global homogenized problem. In order to derive the ho- 
mogenized problem, we need to uncouple problem (|4.4|) . As we can see, (|4.4|) is 
equivalent to the following system f|4.9[) - (|4. 10|) : 
(4-9) 

i,^fr~ dxdt + JJg TXlc aBu ■ dxpidxdtd/3 

b |Du| p ~ 2 Bu ■ dip-tdxdtdP = for all ip 1 G C^(Q T ) O [V At (R t ) ® V di 




xK 



f -OJly(p) J Qt u ■ ip' dxdt + a I {u,(ip Q ,0)) +b I {u,u,(ip o ,0)) 
[ - ' \ =J QT Wl(pf(;U O ))-iP o dxdtfor a \\4, o eC^(0,T)®V. 



It is an easy matter to deal with (|4.9I) . In fact, fix £ G R JVxAr and consider the 
following cell problem: 

(4.11) ^ +J K b\^ + dn(0\ p - 2 ^ + d^))-dwd(3 

I + / ;c a(£ + 0ff(O) -Stod^ = for all w E ^(R r ;^)- 

Assume for a while that the solution to (|4. exists. Then, due to the properties 
of the functions a and b, the linear functional 

«"->• / 6|C + 07?(Of~ 2 (£ + 0w(O)-0«ta/9+ / o(f + fl5r(0) ■ 
«//c Jk 

defined on B p Ap {R T ] B^) is continuous, so that the linear functional to i-> [p^ 51 , to] , 
defined on £div = X>ap(R t ) (g) Vdiv is continuous when endowing £^ lv with the V- 
norm. From the density of £div m ^ P (R T ; B^?) (see Lemma [9]) we get readily 
pf-jr G £>^ P (R T ; (B d ^)'), so that 7r lies in Wdiv Since £div is dense in Wdiv, Eq. 
(14.111) still holds for w G Wdiv But if 7Ti =7Ti(£) and 7r 2 =7r2(£) are two solutions 
of (14. lip then, setting 7T=7ri— 7T2, 

J K Wtc ■ dwdf3 + J K (b(; £ + 9tti) - &(-, £ + <9^ 2 )) ■ dwd/3 = 
for all w G Wdiv 

Taking the particular test function w =ir and using the equality [p^7,tt] = 0, we 
are led to 

f adn ■ dndfi + [ (b(-,£ + dfti) -b(-,£ + dTV 2 )) ■ dndfi = 0, 
Jk Ja(a) 

and using once again the properties of a and &, we get 



[ \dn\ 2 d(3 + v x [ \dn\ p df3 = 0, 
Jk Jk 



which gives dn — 0, or equivalently, D y 7r = 0. It then follows that 7r = since 
it belong to B P AP (R T ; (B^ p Ap (R^)) N ). Returning to the existence issue, Eq. gTT]) 
admits a solution; see e.g., [351 Proposition III.4.1]. 

Now, choosing %p 1 = ip®w in P~9]) with tp G C$°(Q T ) and w G [X>ap(R t ) ® Vdi v ], 
and using the arbitrariness of <p we get the following equation: 
(4.12) 

pm, ^ + / oBu-dtod/3-l- / 6 |Dm| p ^ 2 Du-dwd/3 = for all G £UpQM®Vdr 
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(4.13) 



Coming back to (|4.11[) we choose there £ = Duo(x,t) (for arbitrarily fixed (x,t) £ 
Qt)- Comparing the resulting equation with (|4.12|l and using the density of 
T>ap(J&t) ® Vdiv in £>^ P (R r ; S d - V ) we get by the uniqueness of the solution of (|4.11j) 
that iti = tz(Duq), where tz(Duq) stands for the function (x,t) i-» tt(Duq(x, t)) 
considered as acting from Qt into Wdiv This shows the uniqueness of the solu- 
tion of (14.91) and also allows us to conclude that U\ lies in LP(Qt\ Wdiv) so that 
(iio,Mi) € ^o' P - This concludes the proof of Theorem[SJ 

In order to get the macroscopic homogenized problem, we set, for £ £ M. NxN , 

M(0= [ b\H + dn(0\ p - 2 (l; + dn(0)d(3 

JK 

and 

<= I a(£ + dn(Z))df3. 
Jk 

Next, consider the following problem: Find (1*0,90) solution to 

m v(p)^f ~ div(m£« ) + M(£>uo)) + B(u Q ) + Vq = DJl(pf(;u )) 
div Uq = in Qt 
u = on <9Q x (0, T) 
1*0(2:, 0) = m" (x) in Q. 

Problem (|4.13p is the macroscopic homogenized problem. In contrast to (|4.4[) it only 
involves the macroscopic limit (ito, go) of the sequence (u e , q £ ) of solutions to (|2.1I) . 
Thus it describes the macroscopic behavior of the above-mentioned sequence as the 
heterogeneities go to zero. 

Returning to Problem (|4.13[) we know that < Xft y (p) is a positive constant (see 
Assumption (A5) at the beginning of this section). Moreover one may easily check 
that the functional u h-> 9Jl(pf(-, it)) maps H into L 2 {Qt) n , and is Lipschitz contin- 
uous (this is due to the properties of / and p). Thus, assuming the existence of the 
velocity Uq in the above problem, it comes by [35l Proposition 5] that the pressure 
<7o will exist. We can hence argue as in the proof of Theorem Q] to get the existence 
of a solution to (|4.13l) . We may therefore state the macroscopic homogenization 
result. 

Theorem 7. Assume that (A1)-(A5) hold. Let l + ]frTj < P < °°- For each e > 
let (u e ,q e ) be a solution of (|2.1[) . Then there exists a subsequence of(u £ , q £ ) £ >o (not 
relabeled) which converges strongly in L 2 (Qt) N (with respect to the first component 
u e ) and weakly in LP (Qt) (with respect to q £ ) to the solution of (|4.13[) . Moreover 
any limit point in L 2 (Qt) n X LP (Qt) (in the above sense) of (u e ,q £ ) e> o is a 
solution to Problem (|4. 13|) . 

Proof. Considering (|4.8[) , by a density argument, this equation still holds for 4> in 
L p (0,T; Wq P (Q) n ) x LP(Q T ; W div ). So substituting m = ir(Du a ) in gU and 
choosing the particular $ = (ip , ipi) with xp 1 — and -0 O G C^(Qt) n , we end up 
(using some obvious computations) with the result. □ 
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